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J> ' Abstract 
l> , 

OO , We analyze the interplay between dispersion relations for the spectrum of Hawking quanta 

■ and the statistical mechanics of such a radiation. We first find the general relation between 

the occupation number density and the energy spectrum of Hawking quanta and then study 
several cases in details. We show that both the canonical and the microcanonical picture of 
I the evaporation lead to the same linear dispersion relation for relatively large black holes. We 

also compute the occupation number obtained by instead assuming that the spectrum levels 
out (and eventually falls to zero) for very large momenta and show that the luminosity of black 
holes is not appreciably affected by the modified statistics. 



1 Introduction 



^ ' The gravitational collapse as described by general relativity can lead to the formation of space-times 

■ with peculiar causal structure and profound consequences on the (quantum) matter propagating 

on it. The most striking example is probably that, once the (apparent) horizon has formed, Hawk- 
ing radiation [||| generically sets off Q and the mass of the source should then decrease by this 
quantum mechanical process. The main problem that remains with such a semiclassical picture is 
the determination of the backreacted metric and the corresponding time evolution of the source. 

A different (but possibly related) problem is the role played by very high (trans-Planckian) 
frequencies. In fact, by tracking back in space a photon of frequency u! as measured by a distant 
observer, one immediately finds that its frequency is blue-shifted, in the optical approximation and 
neglecting the backreaction, according to the formula 

^ ~ (^1 - — j ^ ' (1-1) 

where is the radial coordinate at the point of emission. It is clear that uj* is unbounded from 
above if approaches the horizon radius Rh, and one is led to conclude that, in order to have 
UJ finite, co* will soon exceed the Planck mass if ~ Rh as it is expected for Hawking quanta. 
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This is a very strong conclusion, since to study such energetic states one would need a full-fledged 
theory of quantum gravity. 

More recently it was shown that, contrary to the above argument, trans-Planckian frequencies 
do not play a significant role in the Hawking process and the evaporation looks indeed insensitive to 
the presence of a UV (short distance) cut-off (for a review and list of references, see Q). This opens 
up the possibility that the spectrum of emitted quanta be not linear and a non-trivial dispersion 
relation avoid the production of trans-Planckian modes (see, e.g. Refs. ^, |6| for interesting 
proposals). About the origin of the new dispersion relation little is known. One might argue that 
the blue-shift in Eq. (1^) must be corrected for the true (backreacted) metric in the vicinity of the 
(apparent) horizon. In fact, since the Hawking quanta are produced inside the potential barrier 
that surrounds the horizon, a fraction of them gets trapped and forms a (thermal) bath which 
backreacts on the metric. Moreover, the quanta which escape through the barrier do not propagate 
in vacuum since they must cross such a bath 0]. The correct dispersion relation must then account 
for both aspects and, as such, refiects our inability to solve the main problem with black hole 
evaporation. 

In the present paper we shall explore the connection between the dispersion relation and the 
statistical mechanics of Hawking quanta. We shall first work out general expressions in Section |2| 
which we then apply to both the canonical picture and the better sound microcanonical picture in 
Section ^. The starting point of the latter approach is the idea that black holes are (excitations 
of) extended objects (p-branes), a gas of which satisfies the bootstrap condition |^, [l^. This 
yields a picture in which a black hole and the particles it emits are of the same nature and the 
statistical mechanics of the radiation then follows straightforwardly from the area law of black 
hole degeneracy [ pT]] . One obtains an improved law of black hole decay which is consistent with 
unitarity (energy conservation) and no information loss paradox is expected. In fact, black holes 
approximately decay exponentially in this picture, although departures from the canonical behavior 
occur only around (or below) the planck mass ]T0[ |. 

Of course, the statistical mechanical approach is global and does not allow us to fully determine 
the local behaviour of the fields (although some explicit connection with the dynamics of the local 
geometry and the backreaction can be drawn ||l^). In particular, it yields the occupation number 
density of Hawking quanta, but one then needs an extra hypothesis to determine the wave modes 
that lead to such a density. In principle, the new wave modes should describe the propagation in 
the black hole metric with the backreaction included. An hypothesis of this sort was put forward in 
Ref. Q and we shall show that quantitatively negligible corrections to the linear dispersion relation 
predicted by the canonical picture are required by the microcanonical treatment of the Hawking 
radiation, except that there is a natural cut-off at u; = M, where M is the black hole mass. This is 
practically ineffective for the problem of trans-Planckian frequencies since M > 1 (in units of the 
Planck mass, with c = fi = G = 1) for a (classical) black hole in four dimensions. 

In the last part of the paper. Section ^, we shall reverse our line of reasoning and assume a 
dispersion relation in order to determine the corresponding occupation number and compare it with 
the canonical and microcanonical quantities. In particular, we shall choose a spectrum of the type 
proposed in Ref. Q and show that it does not produce appreciable modifications to the luminosity 
of large black holes (in agreement with the general framework of Ref. ) . In Section we conclude 
by mentioning that this result might be significantly modified by the existence of extra dimensions 
1 13 1, as the (microcanonical) luminosity was shown to depend strongly on the dimensionality of 
space-time p 
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2 Occupation number density and wave modes 



An easy and instructive way of obtaining the standard (canonical) occupation number density of 



Hawking quanta is the fohowing |15]. Consider a spherically symmetric four-dimensional metric in 
the Painleve-Gulstrand form 



-c^{r,t)dt^ + [dr - v{r,t) dtf + r"^ dfl"^ , 



(2.1) 



where dQ'^ is the the line element of a unit two-sphere. The metric admits an (apparent) horizon 
if r = Rh exists such that v{Rh) = vh = —ch = —c{Rh)- The surface gravity is then given by 



where 



gH_ 

CH 



1 dic^- t.2 



5H= 2 



dr 



r=RH 



The wave modes 

0(r, t) = A{r, t) exp [ip{r, t)] = A{r, t) exp 



iujt — i / k{r') dr' 



solve the d'Alambertian equation in the eikonal approximation, 

d^Lp d'^(f + ie = , 
provided the wave-number is given by 

k 



UJ 



a {1 + ie) c + V ' 

where a = +1 (—1) for outgoing (ingoing) modes. 

Near the horizon (r ~ Rh), ingoing modes (cr = —1) have wave-number 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



UJ 

'2^ 



(2.7) 



and are defined for r — Rh both positive and negative. Instead, purely outgoing modes (cr = -|-1) 
exist only outside the horizon (r > Rh) with 



UJ 



K{r - Rh) 



(2.8) 



Additionally, in this set of coordinates one can consider "straddling" modes that are defined for 
all values of r > and are swept "downstream" inside the horizon. For such modes one has (again 
for r Rh) 



UJ 



^straddle (?^, t) K, exp iiujt - i - \D.\r - Rh\ 

K 



■K UJ 



e(r - Rh) + exp { +— \ Q{Rh - r) 



(2. 
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where the Boltzmann-hke factor inside the square brackets emerges from the usual analyticity ar- 



gument which relates ingoing and outgoing amplitudes |16] and corresponds to the inverse Hawking 
temperature 



PH — • 

K 



(2.10) 



In fact, the physical vacuum is defined with respect to Astraddle) since freely falling observers should 
not see any peculiarities as they cross the horizon. The Bogolubov coefficients of the transforma- 
tion from the basis {i;^>straddie} to the basis {</>in, 0out} are then simply given by the (normalized) 
amplitude of the ingoing (Mn) and outgoing (A'^out) parts of the "straddling" modes. From Eq. ( ^ ) 
one finds 



|iV. |2 = qI^h^ 



out I 



The (wronskian) normalization condition 

|Mn|'-|iVout|' 

then yields the thermal occupation number density 



(2.11) 



(2.12) 



\N. 



out I 



(2.13) 



for the outgoing Hawking quanta. To summarize, one has obtained the occupation number density 
as a consequence of the near horizon geometry. 

By reversing the above argument, one could in principle assume a specific function for the occu- 
pation number density and then reconstruct the related possible metrics. If the exact n were known, 
one would obtain some insight for the metric which takes the backreaction properly into account. 
It is therefore useful to rewrite some of the above expressions in terms of a (this far) unspecified 
function n{uj). In particular, by replacing with n one obtains new Bogolubov coefficients such 
that 



lAT. |2 = eWi+"-'H] \M 



2 

out I ) 



and the backreacted "straddling" modes are determined as 



"straddle 



(r,t) 



exp {lujt In 

^1 27r 



l + n 



[UJ 



In |r — R}{\ 



e(r - Rh) + VI + ""H^) ®{Rh - r) 



again in the vicinity of the (apparent) horizon (r ~ Rh)- 



(2.14) 



(2.15) 



3 Canonical and microcanonical dispersion relations 



For the following, it is useful to introduce the dimensionless k = {r — Rh) /^out- From Eq. ( |2.15 ) 
and any number density n one finds that, for r > Rh, 



k » In [l + n^^{LL>) 



(3.1) 
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Figure 1: Behavior of um (solid line) and rip (dotted line) for M = 10. Note that um = for 
u > M = 10 while is (formally) defined for all values of w > 0. 

This relation is in general difficult to invert, depending on the form of n(io). For n = np, Eq. ( p.l| ) 
reduces to k Phlo and, for fixed values of r > Rh, one finds 

i « ^ = ^ . (3.2) 

dk (3h SttM ^ ^ 

We now recall that the occupation number density in the microcanonical ensemble is certainly 
a better approximation to the (unknown) exact expression Q. It can be obtained directly from the 
area law without solving for the wave equation and is given by [|8|, |9| 

k exp(4vrM2) ^3 3^ 



w > M 



where M = 1/4 k is the (instantaneous) black hole mass and the function C is a (unknown) factor 
which might account for high energy corrections coming, e.g., from string theory. In the following 
we shall set C ~ 1 unless differently specified. We also note that there is a natural cut-off at uj = M 
(> 1). For a comparison between hm and n^j see Fig. |l]. 

For n = um , k is a rather complicated function of lo. However, one can make some approxima- 
tions on considering that we are particularly interested in the high frequency regime, i.e., uj ^ 1. 
For M ^ 1, it is useful to consider M as a large integer and, for 1 — e < u; < 1 (with < e ^ 1), 
one can easily compute 

duj ^ UM (1 + riM) 
dk uj{dnM/duj) 
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M 


nM 


duj 


di^^ 


1 

SttM 


1 


3.5 X 10"^ 


-AA X 10"^^ 


8 X 10"^ 


4 X 10"^ 


10 


2.0 X 10^^"'' 


-4.5 X 10"^"^ 


4 X 10"^ 


4 X 10"^ 


10^ 


9.0 X IQ-^"^^' 


-2.3 X lO^^'J^^ 


4 X 10-"^ 


4 X lO""* 


10^ 


2.8 X lO-l'^yl" 


-7.0 X 10-™ 


4 X 10"^ 


4 X 10"^ 



Table 1: Occupation number density in the microcanonical picture for lv = 1 and comparison 
between microcanonical and canonical dispersion relations for uj ^ 1 and various values of the 
black hole mass M. 



In this range 



M 

^-4ttM'^ ^ 

1=1 



8tt Mw 



and 



duM 



M 

E 

1=1 



87r(;a;-M)e^^(*^-''^)' ~ -8^Me" 



-Ml 



dnp 
dw 



from which 



duo 



1 



SirM 



(3.5) 



(3.6) 



(3.7) 



in agreement with Eq. (3.2). The above turns out to be a rather good estimate for M ~ 10 and 
greater, as one can check numerically (see table Q). 

For small values of M (~ 1), one must properly take into account the finite sum appearing in 
Eq. (p.3|). The result of a numerical analysis is shown in Fig. ^ for M = 1 and in Fig. ^ for M = 10. 
As was expected from the figures given in Table ||, the dispersion relation for M = 10 is visibly 
linear and does therefore not differ from the canonical picture. For M = 1 the curve departs from 
linearity and turns upward for lj approaching the Planck energy. 



4 Occupation number from dispersion relations 

Upon solving Eq. ( |3.1| ) for n{uj) one obtains 



n = ^, (4.1) 

e'^ — 1 

which is uniquely defined only for intervals of uj where the function k{uj) is single valued. This 
seems to be true for the two cases inspected in the previous Section, and would also hold for a 
spectrum which goes asymptotically constant for large k [Q]. However it does not apply to the 
Epstein functions suggested in Ref. P]. The latter can be considered as an extension to all values 
of k of the spectrum studied in Ref. 

12 / 12 \ 

(4.2) 
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Figure 2: Microcanonical dispersion relation for M = 1. 
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Figure 3: Microcanonical dispersion relation for M = 10. 
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which is defined only for k < A and is presumably meaningful only as the next-to-leading order 
expansion at small k of the correct dispersion relation. 

We shall here consider a particular case of the family of functions studied in Ref. P], namely 



UJ 



+ 



(4 - 2 e) e^/^^c: 



(4.3) 



where kc determines the location of the maximum of oj. We also demand that n{uj) ~ np{ijo) 
for a; <C 1. This uniquely determines the coefficients k^ = (3h and e = from equating the two 
lowest order coefficients in the Taylor expansion of Eq. (4.3) near /c = to the right hand side of 
Eq. (|3.2| ). The constant kc can be fixed by requiring that the maximum of uj is close to 1. Since the 
derivative of uj{k) with respect to k vanishes for k = km ~ (5/2) kc, from uj{km) = 1 one obtains 
kc ^ (3/4) /3h and k^ ^ (15/8) /^h- Finally, 



k , ( 2k 



(4.4) 



In Fig. 1^ we plot a comparison between the dispersion relation (4.4) and the canonical dispersion 
relation ( |3.2D for M = 10 (we recall that the same dispersion relation follows from the microcanon- 
ical picture for such a large mass as shown in Table |l] and Fig. |3|). In Fig. ^ we then display a 
comparison between the corresponding occupation number densities n and n^. We note that in 
the range < w < 1, um > '^/3 while n < n^. One might from this infer that the microcanonical 
number density must be corrected for w ~ 1 by a suitable C(a;) in Eq. (^I^). 

Having determined a novel occupation number density, we can now proceed to estimate the 
corresponding luminosity for an evaporating black hole, which can be formally written as 



C = A 



T{uj) n{uj) uj^ duj , 



(4.5) 



where P ~ 1 is the grey-body factor and A = IGvrM^ the horizon area. We now note that, since 
km ^ 1 for a black hole of mass M > 1, the number density for momenta k > km is highly 
suppressed by the form of Eq. ( |4.1| ). Further, (jj{k) vanishes exponentially for k ^ km- One can 
therefore neglect the contribution of such modes, use the number density in Eq. (|4.1| ) in the range 
< k < km and approximate the luminosity as 



f3' 



H_ 

4 



n{uj) dijj 



4/32 



H JO 



sech 



2k 



2k 

mi 



tanh 



2k 

mi 



k^ dk 



M- 



(4.6) 



where km is the value of k at which the term in square brackets (proportional to dco/dk) vanishes 
and the last line follows from dimensional analysis. The above Eq. ( [4.5| ) is just the standard 
canonical result |||]. The integral can also be estimated more precisely by changing the integration 
variable, 



4 



sech^(x) [1 — X tanh(x) 



dx 



62 



(4.7) 
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Figure 4: Comparison between the dispersion relation (|4.4D (solid line) and the canonical dispersion 



relation (3.2) (dotted line) for M = 10. 
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Figure 5: Behavior of the occupation number density n (solid line) corresponding to the dispersion 



relation (4.4) compared with (dotted line) for M = 10. 
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Figure 6: Black hole luminosity determined by the Epstein dispersion relation ( f4.4D (solid line) 
compared with the canonical luminosity (dotted line). Vertical units are arbitrary. 

where Xm = x{km) ~ 1-2, and performing the integration numerically. The result is displayed for 
comparison with the canonical luminosity in Fig. |^ which shows that Eq. ( [4. 61) is indeed correct to 
an exceedingly good approximation. We finally note that, since the microcanonical luminosity does 
not significantly differ from the canonical expression for four-dimensional black holes (necessarily 
with M > 1 [0), the luminosity ( [4. 6] ) computed in this section is also in agreement with the 
microcanonical result. 

5 Conclusions and outlook 

In this paper we have studied how the statistical mechanics of the black hole evaporation is affected 
by the high energy behavior of Hawking quanta in four dimensions. We have found that one can 
consider large deviations from a linear dispersion relation at near Planckian frequencies without 
changing the luminosity of a black hole. In particular, the fact that the new dispersion relations 
advocated in Ref. Q| do not change the luminosity of a black hole with M > 1 is a direct verification 
of the general framework described in Ref. [^]. We have, however, not attempted at any estimate 
of how such modifications affect the laws of black hole thermodynamics, nor whether they can be 
indeed derived from a fundamental theory. 

Although we have not explicitely considered black holes lighter than the Planck mass in this 
paper, we suspect one would obtain different results for those cases [jlO|. Such objects are outside 
the domain of classical general relativity in four dimensions, since their Compton wavelength would 
be larger then the horizon radius. However, with more than four available dimensions [^] black 
holes could exist with Mq < M < 1 (where Mq is of the order of the fundamental mass scale 
of gravity). Further, the scale below which microcanonical corrections to the luminosity become 
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significant in that scenario is given by the critical mass (much larger then the Planck mass) 



above which the black hole starts to behave hke a purely four-dimensional object |14]. It then 
follows that for black holes with Mq < M < one expects that modification of the dispersion 
relation for frequency 1 < a; < Mc indeed affects the luminosity. We hope to extend our analysis 
along this line in a future publication. 
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